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Abstract 

We present the theoretical approach to study the unconventional 
Josephson junction (uJJ) made by putting the non-superconducting strip 
on the top of superconducting strip. We work in the framework of the 
Ginzburg-Landau, Bogoliubov de Gennes and Usadel equations. We solve 
the non-linear partial differential equations numerically for few simple 
cases. We obtained the eigenenergies of the uJJ by means of combined 
GL and BdGe method for the simplest case. 

We review the similarities and new aspects of uJ J with currently known 
Josephson junctions. We predict the occurence of the physical effect, 
which we call the topological Meissner effect. Basing on the obtained re- 
sults and current knowledge on Josephson junctions we point the future 
perspectives of the research on uJJs. 
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1 Motivation and general overview 



The unconventional Josephson junction (uJ J) can be defined as by placement of 
non-ferromagnetic or ferromagnetic (non-superconducting material) on the top 
of superconducting strip, (as depicted on the Fig.l) The superconductor can 
have order parameter of any symmetry type as s, p, d, f for the case of singlet, 
triplet or mixed case. 

The technological process necessary to produce uJJ is quite simple, since we 
have to evaporate (or bring in another way) the non-superconducting material 
on the top of given superconductor strip. Because of the simplicity of the 
technological process of production of single uJJ we shall be able to produce 
circuits of high or medium level of integration. This has particular importance 
for the d-wave superconductors since some of them can be superconducting in 
the temperature of liquid nitrogen. Such view was suggested by the international 
patent by L.Gomez and A.Maeda [0] that defines the physical structure, which 
we name as uJJ. Richness of physics in such structure is considerable. 

UJJ can be used possibly in the circuits for THZ electronics or superconduct- 
ing qubits. If the ferromagnetic strip is being used, it shall be possible to build 
the tunable Josephson junction with the magnetization of the ferromagnetic 
material as the tunning parameter. Then uJJs should serve as the tunable de- 
tector of electromagnetic radiation or non-linear SQUID. Such SQUIDs should 
be constructed from single or two uJJs and are expected to show highly non- 
linear behavior with the respect to the external magnetic field or with respect 
to the current flowing via the junction or both these factors. Especially inter- 
esting behaviors of non-linear SQUID is for the case when two uJJs strongly 
couple via the ferromagnetic stripes. In such case for the certain regime SQUID 
physical structure should be preventing entrance of external magnetic field to 
the interior of SQUID for certain values of this field. 

Also nonlinear uJJ SQUIDs built from one or two unconventional Josephson 
junction should have higher sensitivity to the external magnetic field. Par- 
ticularly under certain conditions the unconventional Josephson junction shall 
be capable of detection fraction of quantum fluxon. Full theoretical confirma- 
tion of the occurrence of such expected phenomena requires conduction of 3 
dimensional numerical simulations of Ginzburg-Landau equation for the certain 
configurations of the external magnetic field. This task is beyond the scope of 
this paper and will be the subject of the future work. 

In this paper we use mainly the approach based on the combined Ginzburg- 
Landau (GL) and Bogoliubov de Gennes (BdGe) formalisms to get the proper- 
ties of the uJJ. Such choice is due to the fact that relaxation method allows for 
the quick solution of the Ginzburg-Landau equation in quite many cases what 
was confirmed in the conducted numerical simulations. 

On another hand, the numerical solutions of the BdGe equations are more 
demanding, since they involve the necessity to find n eigenenergies and n com- 
plex eigenfucntions what is quite difficult especially in the 2 and 3 dimensional 
case. 

Other formalisms are also capable of even more detailed description of the 
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unconventional Joscphson junction. However the prize to be paid for the more 
detailed answers is the significant increase of the complexity of differential equa- 
tions that have to be solved. Therefore application of more advanced formalisms 
than GL and BdGe remains the subject of future studies already being under- 
taken. 

2 Mathematical statement of the problem 

In order to approach the complex case as the solution of the BdGe and GL 
equations for the uJJ we start from the simpler cases which also occurs in 
the system. Let us consider the case of the rectangular shape ferromagnetic 
or non-ferromagnetic non-superconducting bar on the top of superconductor. 
As it is well known the superconductor has special type of physical proper- 
ties compared to the non-superconducting (normal) state of the same material. 
However the proximity effects modifies properties of superconductor (sc) and 
non-superconductor (nsc) material. 

In general situation dealing with the ferromagnetic material on the top of 
superconductor we have the ferromagnetic material to be in the one among 
many possible magnetized states what brings many mathematical cases to be 
accoimtcd. 

The simplest starting case is the whole system in the normal state. Let us 
assume that the ferromagnetic material will have the similar conductivity as the 
superconducting material in the normal state. Then the current flowing in the 
region underneath the ferromagnetic bar changes the magnetization of the fer- 
romagnetic bar, which becomes nonuniform. The magnetization of ferromagnet 
can be assumed to be uniform only for the case of small thickness of ferromag- 
net. The magnetic field in the ferromagnetic bar exerts the Lorentz force on the 
electrons and holes moving in the normal state. Electric current flowing via uJJ 
generates own magnetic field acting on current carriers by itself and acting on 
the ferromagnet by changing its magnetization. By certain magnetization ferro- 
magnetic acts back on the moving electrons and holes. Therefore the resistance 
of the uJJ structure should depend nonlincarly on the magnitude of the electric 
current, magnetization of ferromagnetic bar in the absence of external electric 
current flowing via the system and on the presence of the external magnetic 
field. 

Therefore for the given geometry of the normal and ferromagnetic material 
the determination of current- volt age characteristics is highly non-trivial and 
cases of ac and dc current or voltage shall be accounted. 

Additional complication is the possible anisotropy of the resistance in the 
given structure, what is the case of the d-wave superconductors and is also due 
to the evaporation of the non-superconducting material on the top of supercon- 
ductor. 

The situation is less complicated for the case of isotropic s-wave supercon- 
ductor. 

D-wave superconductor is a kind of ceramic material in the normal state. 
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Therefore if the uJJ is in the normal state the current will tend to flow via 
the ferromagnetic material, which is less resistive than the superconductor in 
the normal state. It is expected to be for the case of ferromagnetic and non- 
ferromagnetic non-superconducting material on the top of superconductor. 

2.1 Case of rectangular shape superconductor of infinite 
height with dimensions a,b in GL approach and non- 
zero current flowing via the system for s-wave super- 
conductor 

In order to approach the consideration of all cases for the uJJ with s or d 
wave symmetry of the order parameter, it is necessary to consider the simpler 
situation (limiting case), when we describe the state of superconductor with 
total current flowing via as I, when it is very far from the non-superconducting 
strip. 

Therefore we study rectangular shape superconductor of infinite height and 
a and b dimension in its smallest cross section, with total ciuTcnt flowing via the 
system equal as I. We look for the distribution of the order parameter, vector 
potential and current density in the structure. The simplest way to approach 
this problem is to state and solve GL equations, which are only valid for the 
case of temperature of uJJ close to critical temperature Tc- Having the solutions 
of the GL equations, we might look for the more detailed physical properties 
as by solving BdGe, Usadel, Eilenberger or Gorkov equations for this system. 
However those formalisms are significantly more complex and more difficult to 
be solved, so initially we will concentrate only on GL approach. 

In the given system we use Cartesian coordinates as the most suitable for 
its description and set the z axis to be in the direction of the height of the 
rectangular shape superconductor (sc) of s-wave symmetry. Since height is 
infinite the system is only invariant under shift in the z direction. The system 
physical properties are therefore parameterized only by x and y coordinates. 
Because of translational symmetry we expect the electric current to be flowing 
only in z direction. Using London limit we have current density proportional 
to vector potential. Therefore we have only z nonzero component of vector 
potential. The boundaries of the rectangular sc are 

A = {xA,yA) = {x = a/2, y e (-6/2, b/2)),B = {xb,Vb) = {x = -a/2, y e (-6/2, 6/2)) 
C = {xc, yc) = (2/ = 6/2, X G (-a/2, a/2)), D = {xD^yo) = {y = -6/2, x G (-a/2, a/2)) 
We have the complex scalar fleld describing the order parameter given as 
ip{x,y,z) = \ip{x,y)\e'^cA''^^'y)i^-^o) 
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where zq is fixed value. At the boundaries superconductor-vacuum for the su- 
perconducting rectangular suspended in vacuum we have 

Av.(^) = o,|:«B) = o,A^(C) = o, 1^(0 = 

Intuition tell us that the order parameter in the middle of the superconduct- 
ing rectangular {xo,yo) has the maximum value and should be decreasing when 
we arrive to the boundary. 

Thus we have 



^^\'^{^^y)x=xo,y=vo \ ~ 0' ^y\'^{^'y)x=xo,v=yo \ ~ ^ 

GL equation can be written as the following 



aV + /3|V'l'V'+^(PcV') = 

We assume that a has constant value in the given rectangular and zero 
outside. We have the canonical momentum operator pc given as Px,c = ~ 
^Ax) and = p^ ,, + Py ,, + p^ ,,. We set gauge to be VA = and obtain the 
electric current density in GL formalism as 

j = Vx(VxA) = V{VA)-V'^A = -V'^A = ^{tp^Vil>-iljVil>^)-cA\^\'^ = csAlipf 

where C3 is constant depending on the universal physical constants. Because of 
non-zero z component of A current density is j{x, y, z) = jz{x, y). Therefore we 
have the set of 2 dimensional non-linear partial differential equations for the real 
value functions \ijj{x,y)\ and Az{x,y), where \^{x,y)\ is always non-negative. 

-(^ + ^)^^(^'y) = ciAz{x,y)\il){x,y)\'^ 



^{x,y,z){a+p\i;{x,y)\^)-^{^+-^)^{x,y,z)+^{~~A,{x,y 
The last equation is as 

!P (f (P 

ipix, y, z){a + P\tp{x, y)f) - + ^)'^(^' z) + C3A;,{x, yfipix, y,z) = 

where C3 is constant. We have 
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The great complexity of equations can be simplified if we set ^ = 2:0 so we 
obtain 

We might attempt to solve two equations simultaneously or try to reduce 

two equations to the one more complicated equation. 

We observe that we can write the magnitude of the order parameter only in 
the terms of vector potential Az{x, y) what is given below. 



\i^{x,y)\ = 



ciA^{x,y) 



Thus we have the complicated two dimensional equation for the single real 
value function(the z-th component of vector potential that has certain con- 
straints) of two coordinates as . 



To solve such highly non-linear equation analytically or numerically is uneasy 
task. 

Additional constrain for the system is the value of the total current flowing 
via the system as 

r— r- 

dxdyj^{x,y) = I 

where I is the given value of total electric current flowing via the system. 

Experimentally we can set certain value of total electric current I, which is 
below Ic{T). The dependence of order parameter is accounted by the value of 
a{T). It should be underlined that in the superconducting state a is negative 
while 13 is positive. In the normal state both a and /3 are positive. 

The solutions of the stated two dimensional GL equations is not straightfor- 
ward. To exercise our method we will flrst solve the more simplified case and 
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then solve the case currently described. The simpler case under the considera- 
tion will be the reduction of dimensionality in the coordinates. If we consider 
indefinitely long superconducting cylinder of given radius R, with certain cur- 
rent flowing via the system that we have to describe the system by only one 
coordinate r. 

2.2 The solutions of the GL equation for the infinitely long 
superconducting cylinder of radius R with nonzero 
current flowing via the system 

The system has two symmetries: translational in the z direction and rotational. 
This simplifies the GL equations so they become dependent only on r. 

Similarly as before we have the electric current flowing in the z direction so 
jzir) electric current density component is nonzero and hence ^z(r) is nonzero. 

We expect the order parameter in the system to be of the following form 

where zq is the any fixed value of z coordinate. 

One of the additional constrains is that the total magnitude of the electric 
current flowing via the system is given as 

1= rdrjzir) 
Jo 

In cylindrical coordinates we have given the boundary condition as 

-J^|V'(r)|.=fl = 
ar 

Intuition tell us that the order parameter in the middle of the superconduct- 
ing cylinder has the maximum value and should be decreasing when we arrive 
to the boundary. 

Thus we have 

|;|^(r)|,=o = 

We obtain GL equation as 

aij{r,z)+l3ij{r,z)\i;{r)\''-^{-^+i-^))iPir,z) + ^i-^^ 

2m r ar dr'^ 2m i dz c 

Similarly as before a is constant inside the cylinder and is zero outside. 
This equation after taking derivatives and simpliflcation is of the form 
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+ i{iPir,z)-^A,{r)^{z- Zo)+ 
Therefore finally we obtain 

Id d^ 

-(^A.(r)2(^ - ^o)')V(r-) + iCV-Cr-, z)^A,{r)^{z - zo)+ 

-^7/.(r,z)A4r)2ci =0 

This equation can be separated to 2 equations on two real-valued functions. 
Particularly if z = zq we obtain the equation 

a|V(r)| + mr)\' ' + (^))IV'(r-)| = 

Its solutions shall be valid for any z since system has the translational sym- 
metry in z direction. 

It turns out that one of the most effective method to solve this equation is 
the relaxation method. 



2.3 Special case of the vortex presence in the supercon- 
ducting cyhnder of the infinite height and certain ra- 
dius r 

We might assume the existence of the vortex exactly being placed inside the 
infinite-length superconductor so the translational and rotational symmetries 
are being preserved. In such situation we have nonzero A^{r) and Az{r). 

Then there occurs two minima of the magnitude of the order parameter at 
r = and r = R. We have two non-zero currents components in the system 
as component j0(r) supporting the existence of vortex in the sc cylinder and 
component jzir) responsible for the observed effective current I in the system. 
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2.4 General scheme of the relaxation method in the solu- 
tion of GL equations 

The presented cases of the GL equations can be solved numerically using very 
simple rule, as by the minimization of the free energy functional, which is the 
basic assumption in the derivation of the GL equations. The final solution, 
which is the configuration of the fields {tp, A) (effectivly 5 scalar fields) in the 
given space, which can be 1, 2 or 3 dimensional with certain boundary condition 
fulfills the equations. 

—F[ilj[x,y,z],il;[x,y,z]\A[x,y,z]] = 

and 

—F[tp[x,y,z],tp[x,y,z]\A[x,y,z]\ = 

To approach the solution wc need to make the initial guess of the field 
configuration on the given set that corresponds to the certain physical intuition. 
The initial guess should be not so far from the final solution. 

Having the initial guess we perform the calculation of fields changes on the 
given lattice 

S S 
—F[il)[x,y,z],il)[x,y,z]\A[x,y,z]\ = -—il)[x,y,z] 

s s 

—F[ip[x,y,z],ip[x,y,z]^,A[x,y,z]] = -—A[x,y,z] 

In the simulation cases St — At is usually fixed and shall be not too big since 
it might affect the numerical stability of the algorithm and not too small since 
wc would like to conduct the numerical GL solution in finite time with finite 
accuracy. To see concrete numerical examples please refer to figures 7,8 and 3 
and 4. 



2.5 Infinite lenght rectangulcir shape superconductor of 
dimensions a,b in terms of BdGe 

Wc assume the crystal isotropy(what implies the assumption of the lack of 
dependence of the effective electron and hole mass on direction) and s-wave 
order parameter and have the Hartree potential of the form as given below, 
where Vq accounts for the energy necessary for the electron to leave the crystal. 

Vix, y) = Vo8{x - a/2) + Vo8{x + a/2) + Vo5{y - h/2) + Vo5{y + 6/2) 
The Hamiltonian of the free particle is then given as 
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The normalization condition for Un{x,y) and Vn{x,y) functions gives the addi- 
tional constrain as 

/•a/2 rb/2 

/ / dxdy\vn{x,y)\'^ + \un{x,y)\'^ = 1 

J -a/2 J-b/2 

Hu„{x, y, z) + A{x, y, z)vn{x, y, z) e„u„(a;, y, z) 

-H^Vn{x, y, z) + A(a;, y, z)^Un{x, y, z) = e„u„(a;, y, z) 

Wc set the effective coupling constant to be constant inside the supercon- 
ducting crystal as Vx and outside. We set the Fermi-Dirac distribution function 
to be /(e„) = l/(+e^''/'^'^) so we have 

l^{x,y,z) = -Vi^Un{x,y,z)Vn{x,y,z)^{l - 2f{en)) = A{x,y)A{z) = 

n 

n 

We work in the gauge so V A = 0. 

We have given the additional equation for vector potential as 
Mz{x,y) = V'^A;,{x,y) = ^(|u„(a;,y)|^ - |v„(a;, j/)p)(A;^,„ - (2e/c)A^(a;, j/)) 

n 

The last equation can be written in terms of the matrix and eigenvalue 
equation. Here the eigenfunctions are A. 

The last constrain to be fulfilled is that total current flowing via the system 
is I and given as 

/a/2 ph/2 
/ j{x,y)dxdy = I 
-a/2 J-b/2 

2.6 Infinite length superconductor beir of a,b dimensions 
in terms of Usadel equation 

The Usadel formalism is capable of dealing with the situations when the su- 
perconductor is dirty. Here we assume that there is no current flowing via the 
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system. Wc assume the following parametrization function 0„(a;, j/) of normal 
G and anomalous F propagators given as 

G(x,y)n = cos{en{x,y)),F{x,y)n = sin{Qn{x,y)) 
We can also use Riccati parametrization given as 

Gn{x,y) = — J r — -, —,Fn{x,y) = — ^ — 

1 + an(x, yjttnix, y)* 1 + a„(a;, y)a„{x, y)* 

where an{x,y) is the complex function. 

Both parameterizations automatically fulfills the normalization condition for 
normal and anomalous propagators given as 

Fr,{x,y)Fr,{x,y)^ + \Gn{x,y)f = 1 

Here a;„ = ■jTT{2n + 1) is the Matsubara frequency and n is the integer 
number. 

The boundary conditions are as following 

Ifm = 0, iKiB) = 0, = 0, |f„(D) = 

We have the Usadel equation as given below 



2A{x,y)G„{x,y) = -{Fr,{x,y)V^G„ - Gr^V^F^,) + 2uMx,y) 

This equation is further simplification of the Eilenberger formalism, which is 
the specific case of the Gorkov formalism. Here D is the constant that is param- 
eterized by the relaxation time of the scattering of electrons on the impurities. 

We have also the expression for the density of states given as 

Nix,E„) = N{0)Re{cos{en{x))) 

where the -A''(0) is the DOS at the Fermi level. 

If Fn functions are small as in the case of the superconductor close to the 
critical temperature Tc and G„ functions are close to 1. Then simplified and 
linearized Usadel equation is of the form 

2A{x,y) = -£.y^Fn{x,y) + 2w„F„(a;,y) 

The variation of A hmction compared to F variation it is small. Therefore 
locally A can be considered to be constant. We introduce here the function 
Fi{x,y) = 2u>nFn{x,y) — A{x,y) Such equation is the easy to be solved. For 
example it can be written in the matrix form and then look for the eigenvalues 
and eigenfunctions of this matrix. We have the self-consistency relation in the 
way as 
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+ 00 

A{x,y)='KT ^ F{x,y,u„)n 

n— — oo 

We need first to have the order parameter distribution as given initially 
by GL formalism. Then we obtain n F„ functions. Having those functions 
we obtain the new order parameter. We continue this process until the order 
parameter will not change in the next iteration. This means that self-consistency 
is obtained. In such way we obtain the corrections of the order parameter to 
the GL equations. 

If we assume that there is non-zero current existing in the system in z direc- 
tion then it modifies our parametrization and the equation of motion as by 



G{x,y)r, = cos{Q^{x,y)),F{x,y,z)^ = sin{Qr,{x,y))e'^cA.(«:,y)(^-^o) 

2.7 Simple method of solving GL, BdGe, Usadel equations 
on sc square 

We present the numerical method of solving nonlinear differential equation of 
the form 

da^ + = Sifi^^, y)) 

where g{,) is a function of function f(x,y). This method can easily be gen- 
eralized to 3 dimensional case (equation of the cubid). 

We assume that the superconductor is a square what imposes certain the 
symmetry conditions on the GL solutions. We will use 2 dimensional grid rep- 
resented by the table where the x direction is the table vertical direction and y 
direction is the table horizontal direction. 

We use the simple second derivative approximation and set Ax = Ay = Ah 

as 



. , , / X ip{xi+i,yj) - 2'tp{xi,yj) +tp{xi_i,yj) ^ ip{xi,yj+i) - 2ip{xi,yj) + ip{xi,yj_i] 
^dx^^dy^'^''^'^''^''^ Ax^ + Ay^ 

= -^{-4:i'{xi,yj) + 'tp{xi+i,yj) yj) ^'(a^i, %+i) + tp{xi,yj-i)) 

We assume that in the geometrical center of the sc square '4>{xo^yo) = a and 
that in the neighborhood places the GL takes values b and c as depicted below. 



c 


b 


c 
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a 


b 


c 


b 


c 
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Having the knowledge of a immediately we get the knowledge of b and c. 
Extending our neighborhood of the initial center point we get the following grid 
of values. 







f 


c 


(>■ 


e 


c 


b 


c 


e 


f 


b 


a 


b 


f 


c 


c 


b 


c 


e 


g 


6 


f 


e 


g 



First we establish values f. Then we compute values of e. After this we 
compute the value g. Having these values computed we extend our initial grid 
in similar fashion and compute all the next values. 

In such way we might determine all values of the GL on the square. The 
last thing we need to check is to confirm that the boundary conditions at the 
sc-vacuum interface occurs. If this is not the case we have to take another value 
in the square center. 

The simple algorithm is also capable of solving more complex systems hav- 
ing certain symmetries, for example the superconducting strip surrounded by 
ferromagnetic layer or non-ferromagnetic layer as given below. 
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We assume that in the presented structures the magnetization of the fe is due 
to the magnetic field produced by the current flowing via the superconductor. 

2.8 Infinite length cylindrical shape superconductor of ra- 
dius R in terms of Usadel formalism 

We assume the following parametrization function 6n(r) of normal G and 
anomalous F propagators given as 

G{r,Un)n = cos{en{r)),F{r)n = sm(e„(r)) 

Such parametrization automatically fulfills the normalization condition for 

normal and anomalous propagators G and F. 
We have the Usadel equation given as 

2A(r)G„(r) = -|(G„(r)(^+ A)i.„(,)_j.„(,)(^+_l)G„(.))+2u;„F„(r) 

After using the parametrization for the case of zero electric current and zero 
magnetic fields in the system we have 
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2A(r)cos(e„)(r) = -|(cos(e„)(r)(^ + A.)sin{en{r)) 

cP d 

-sm(e„(r))(^ + —)cos{@n{r)) + 2w„sm(e„(r)) 

Such equation is second order non-linear ordinary differential equation (ODE) 
and can be solved numerically quite easily. We have the boundary conditions 
given as -^Fn^ii=r{r) = 0. By choosing certain discrete step Ar and by guessing 
the first value at the r=R we move towards the center computing all values on 
the lattice. After reaching the center we compute the derivative. If it is non-zero 
we start guessing from another number until this condition is satisfied. 

The final solution should be of the form that A(r) is maximal at r=0 and 
in monotic way decreases as we move towards r=R. 

Having evaluated @n{r) for r G (0, R) for given n we change index n. Then 
we compute the new order parameter. After some iterations self-consistency is 
achieved. 

We can also account for the electric current density flow, as for example in 
z direction. Then we have 

F{r,z) = \Fir)\F{z) 

where F(z) is the same phase factor as for the order parameter. 

The electric current density flowing via the superconductor is given by the 
formula 

j,{r)=nTDN{0)2iey'{{F^{r,z)\~--A,{r))F^{r,z))-{F^{r,z){- 

^ — ^ I az c I az c 

n 

= -7rTD^c?|F(r)„pA,(r) 

n 

where; N(0) is the density of states at the Fermi level and Ci is the constant 
depending on the fundamental physical constants. 

To flnd solution we need to deal with two coupled differential equations. 

2.9 Infinite lengtli cylindrical shape superconductor of ra- 
dius R in terms of BdGe 

We have the Hartree-Fock potential given as 

V{r) = Vo6{r - R) 
We have the order parameter as given by 

A{r,z) = A(r)A(z) = | A(r)|e*2t^'('')(^-^'') 
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M„(r,2) = |u(r)„|e^t^4'-)(^-^o) 

We have the BdGe equations in terms of 

(f d 

^((^ + ^) + {kz{r) - {e/c)A^{r)f)un{r, z) + A(r, z)vn{r) = CnUnir, z) 

2^^^d^ + + ('^z (0 + {e/c)A,{r)f)vn{r, z) + A(r, z)^Mn(r) = e„i;„(r, z) 
The self-consistency equation gives 

|A(r)| = -Vr J2 \un{r)\K{rm - 2/(e„)) 

n 

We have the Maxwell equation correlating the rotation of the magnetic field 
with the current density as 



H^\un{r)\))e''^^'(^^^'-^°\^{z-zO)j^A,{r)+{\^^^^^ 

To get the simplification we set z = zq what significantly simplifies our 
equations. 
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2.10 Issue of dimensionality of uJJ 

The investigation of the uJJ shaU be overtaken basing on the siniphfication 
of 3 dimensional problem to 2 dimensional problem and assumption that the 
lenght of non-superconducting strip, which is placed on the top of infinite area 
superconductor of finite thickness. In such case the relaxation algorithm was 
applied to solve the GL equation in the absence of magnetic field. After we have 
obtained the order parameter, we plug it to the BdGe equations and look for the 
eigcnencrgies of the system. It should be underlined that the reduction of the 
3D dimensional structure to 1 dimensional problem seems to be not appropriate 
since in such case we loose many physical properties of the system. 

2.11 The simplest case of uJJ in GL picture solved by the 
relaxation algorithm 

We simply the problem of uJJ to two dimensions and consider the non- 
superconducting strip to be non-ferromagnetic. If there are no currents in the 
system, the GL equation is significantly simplified and il){x, y)is real and is given 
as. 

aip{x,y) + pi^{x,yf - _(_ + -j^mx,y) = 
We parameterize the system by a,b,c,d numbers. 

a = const for x G (—a, a), y G (—6, b), E=(-a,-b), F=(-a,b), C=(a,b), D=(a,- 
b), A=(-c,b), B=(c,b),Ai=(-c,b+d),Bi=(c,b+d). 

We have the following boundary conditions in the system 

^i^iFE) = 0, £^iCD) = 0, ^mD) = ^^i^FA) - Q^^BC) = 0, 
^^/KAB) = j^ij{AB), i^iMBi) = Q,±i,(AAi) = 0,£i^{BB,) = Here h 
is real constant value that describes the behavior of order parameter at the NS 
interface. 

For some cases we can assume that the thickness d is infinite. We can also 
assume for some cases that a — )■ ±00. 

The Fig.l3 presents the geometrical parametrization of uJJ. 

2.12 The 2D case of uJJ in BdGe picture 

We are given the order parameter in the form as 

A{x,y) = |A(a;, ?/)|e'^^''^^'^'^^~^°^+'^^''^^'^'^^~^°^ 



Un{x,y) = \un{x,y)\e 



iA^{x,v){x-xo)+iAy(x,y)(v-vo) 



Vn{x,y) = |u„(a;,y)|e'^»'(^'^)(^-^°)+*^''(^'^'(^-^''' 

where xq and yo are fixed values of x and y coordinate. We have the Hamil- 
tonian of the form 
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2m dx'^ dy'^ 

where H is the Hartree-Fock effective potential for electron or hole. 
We define state vector as 




and we have 

- 0, £.P{CD) = 0, ^^i^iED) = A) = 0,^^i;{BC) = 0, 

^^(AiBi) = 0,£mAi) = 0,£mBi) = 

and u'^{SN+) - <(5iV+) = Qu{SN) and v'^{SN+) - v'^{SN+) = Qv{SN) 
where SN denotes the interface between Sc-N and Q is constant depending 

on Sc and N. 



3 Properties of tt Josephson junction. 

3.1 TT Josephson junction in the Eilenberger formahsm 

Using the Born approximation we might simphfy Gorkov equations to the Eilen- 
beger equation. We have the following normal and anormal equasilclassical 
Green functions / = /(a;„,wj),g = g{ujn,Vf). Wc have then the impurities 
present in the system on which the quasiparticles are being scattered with mean 
relation time given as r. Let us consider the system as depicted on the figure 
19 where we assume the problem to be one dimensional. We have the following 
equations in Fe region : 

^/x^5=(2T)-i(/</t>-/t</>) 

^f^^f = -(2t^n + ^h)f + iT)~\f < 9 > -9 < f >) 

-Vf^^-^.f' = -(2w„ + ih).f + {t)-\P <g> -g < f >) 
Here brackets denontc the averaging over the Fermi surface. 

j = 27rr7V(0)^ < {Vflm{g{ujr,))) > 

n 

Here Vf means the Fermi velocity,-/V(0) is the density of states (DOS) and h is 
the magnetic field strength. 

In Sc region we have the equations as 

dx 

Vfx-^f = -2L0nf + 2Ag 
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Wc use the approximation as hr >> 1 what makes < / >=< >= 0. 
In general case the current-phase relation can be expected to be of the form 
= T,t=ijriSin{n^), where jn is the n-th harmonic of Josephson current 
and ^ is the phase difference between superconductors. 

3.2 TT Josephson junction solved by the means of the re- 
IcLxation method 

We obtain the tt- Josephson junction solved by means of the relaxation algorithm. 

We have assumed the constant magnetization in the system in the direction z. 
The problem becomes one dimensional (we parameterize system by z coordinate) 
and is being solved by the relaxation method. Please refer to the Fig.ll. 

3.3 Ferromagnetic grain on the top of superconductor and 
modification of the relaxation method 

The case of the rectangular shape superconductor of dimensions a,b with non- 
zero current flowing via it has been solved. However in our system, we deal 
with the ferromagnetic material on the top of superconductor. In such case the 
current flowing via the superconductor magnetizes the ferromagnetic element 
on the top of it. Reversely the magnetized material creates additional magnetic 
fleld that modifies the distribution of the current in the system. 

3.4 Ferromagnetic-superconductor structure and further 
modification of order parameter 

The situation becomes more complicated when we have to consider not only one 

ferromagnetic grain but entire assembly of ferromagnetic grains. Then we have 
to take into account not only the interaction of each grain with superconductor, 
but also the interaction of grains with themselves. 

3.5 Resemblance of the unconventional Josephson junc- 
tion with existing systems. 

It should be pointed that very essential for the Josephson effects are the dimen- 
sions of the system and two characteristics numbers: superconductor coherence 
length ^ and magnetic field penetration depth A. 

If we use the s-wave superconductor those last quantities ^ and A do not 
depend on the direction. However, in the case of the d-wave superconductors 
they depend on the direction and we can distinguish ab direction and c direction 
so we have (ah,Cc and Xah,K- 

We have two main regimes when the size of the non-superconducting element 
is comparable with ^ab or and when it is much bigger. In the first case we 
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will have the signle Josephson junction and in the second case we will have the 
double Josephson junction. The situation becomes even more complicated when 
we use the ferromagnetic material as the non-superconducting material. 

When the uJJ is being placed in the pc^rpc^ndicular magnetic field than it 
shall have certain features of the tt— Josephson junction. 

However if the material is not being magnetized than it has more resemblance 
to a weak link and wc have the system as Sci — Sc2 — 5*01 . 

All type of the Josephson junctions occurring in our system are field induced 
Josephson junctions. In some sense they belong to the new class of Josephson 
junctions. However they can be approximated by already known Josephson 
junctions. We have at least 5 limiting cases given below. 



Case 


Non-sc 


Magnet iz. direct. 


Ext. mag field 


a/b ratio 


Final effect 


1 


Fe 


_L to surf 


_L to surf 


« 1 


ScFcSc 


2 


Fe 




_L to surf 


« 1 


ScFeSc 


3 


Fe 


to surf 




« 1 


SciFe-^Sc2FeiSci 


4 


Non-Fc 






« 1 


Sc-N-Sc 


5 


Non-Fc 




_L to surf 


« 1 


Sc-N-Sc 



It should be underlined that the uJJ shall be approximated by the networks 
of the Josephson junction. The most simple approximation is about one dimen- 
sional network of the Josephson junctions; the more advanced approximation is 
two or three dimensional network of Josephson junctions. 

Some of the Josephson junctions in this network are being coupled. One of 
the simples model representing the coupled Josephson junctions is give as below 



vac 


vac 


^ fe ^ 


vac 


vac 


sc 


set 


sc 


sc I 


sc 


vac 


vac 


<— 


vac 


vac 





x JJt 




X JJ 4- 





It should be underlined that the fringe field coming ferromagnetic bar is 
responsible for creating two Josephson junctions. Such results were considered 
by the first time by by W.Clinton. It shall be underlined that higher external 
currents flowing in x direction via the system should uncouple two existing 
Josephson junctions so there will be no more polarization as depicted above. 



vac 


vac 


fe B O fe 


vac 


vac 


sc current flow direction sc 


sc 


sc B (g) sc 


sc 


sc current flow direction sc 


vac 


vac 


vac 


vac 


vac 



It should be underlined that if the superconductor is the d-wave supercon- 
ductor and if the c axis is in vertical direction, than the last configuration of 
magnetic field should induce pinned Josephson vortices in the superconductor 
area underneath of the ferromagnetic strip. 



21 



3.6 Andreev reflection in the system and Andreev bound 

states 

Because of defect of the order parameter underneath of the non-superconducting 
system we shall observe the occurrence of the Andreev reflection in uJJs. Even 
more complicated behavior shall occur when the non-superconducting strip is 
being magnetized. Then there suppose to occur the currents in the supercon- 
ductor that will try to compensate the occurrence of the external magnetic field. 
The presence of the magnetic field should double the existing quantized level in 
the Si ~ S2 — Si system, where ^2 superconductor has smaller order parame- 
ter that Si singlet superconductor and has the occurrence of fractional triplet 
phase. Andreev reflections are well described in the framework of the BdGe 
theory. In order to determine the existence of the Andreev states in the system 
we have designed new algorithm for solving 2 dimensional BdGe after plugging 
the initial order parameter distribution from GL relaxation method. 

In Sc(L)-N-Sc(R) system in N region there occurs the Andreev reflection, 
which can be described by the equation 

/ Pelectron{E)dl - Phole{E)dl + (j) - l3{E) = 2^71 
JL JR 

where n is the integer number, j3 is the additional phase shift depending on 
the shape of scattering potential A(r), (j) is the phase difference between su- 
perconductors L and R. We can generalize this formula quite easily for 2 or 
3 dimensional case. We can have the spin-flip in Andreev-bound states when 
the magnetic field in the N layer (Fe layer) is not uniform. Such process is 
schematically being described at the figure 14. 

3.7 Algorithm of solving the BdGe equation in 2 or 3 di- 
mensions 

We map 2 dimensional order parameter distribution to the one dimensional vec- 
tor. Basing on this we build 2 dimensional matrix representing discredited 2 
linear BdGe equations among existing 3. Having the matrix given, we determine 
its eigenfuctions and eigen energies. Having the eigen energies and eigenfunc- 
tions as the one dimensional vector, we map it into 2 dimensional space. Then 
by self-consistent equation we obtain the new order parameter in the system. 
We repeat the whole process until the order parameter will not change in the 
next iteration. 

Let us review of the simple 5x5 grid and look for the solution of the BdGe 
on it. 

We have map two dimensional matrix 
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^5,1 


-^5,2 




-^^5,4 




Xi.l 


X4.2 


^4.3 


A'4.4 


X4.5 


X'.A 


A:!,2 


A3,:! 


A3.I 


A,,5 




^2,2 


X2,3 


X2A 


X2,5 




Xi^2 


-^^1,3 


Xl^4 





into one dimensional vector as given below. Then the paramctrization by two 
indices is changed into parametrization by one index. Mapping 2 dimensional 



Table 1: Two dimensional grid for X(^x,y) functions that can represent 
Un{x,y),Vn{x,y) functions 



Xo ^ 


Xi ^ 


X2 ^ 


X3i 


iXr^ 


Xe ^ 


X5 ^ 


Xi 


Xs ^ 


Xg ^ 


Xio — > 


Xii i 


Xu ^ 


Xi4 ^ 


Xis ^ 


X\2 



matrix to one dimensional vector allow us to write 2 among 3 BdGe equations 
in the compact way. Now we need to make one longer vector representing u and 
v functions as 



We need to find now the structure of the matrix A. On its diagonal is has the 

order parameter mapped to the ono-dimensional vector in the same fashion as 
u and V functions. Now we need to represent the action of 



n,i,j-l) 



2m ^ {Un,i+l,j - '2u„^ij + Un,i-lj) {Vn,i,j+1 - 2l^n 

We find the structure of matrix numerically and than we find its eigevalues and 
eigenfunctions. 



3.8 The occurrence of the vortices in the uJJs 

We should account for the occurrence of vortices in uJJs. This is especially to be 
expected when we have the occurrence of ferromagnetic strip on the top of the 
superconductor and the presence of the external magnetic field. The situation 
becomes even more complicated when we consider the d-wave superconductor 
with ab-plane in parallel to the external ferromagnetic strip. Than we will have 
the possible occurrence of the Josephson and Abrikosov vortices and pancake 
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vortices as well. What is interesting, if we switch external magnetic field per- 
pendicular to the superconductor surface that we might have the switching of 
the Josephson to Abrikosov vortices. Since Abrikosov vortices have smaller core 
than Josephson vortices, for certain magnetic field magnitude, they will allow 
for higher superconducting current flow and thus contribute to possible increase 
of the critical current. 

3.9 Vortices in the superconductor-ferromagnet system 
3.9.1 Abrikosov vortex in the superconductor 

We describe the radial defect of the order parameter caused by magnetic field 
punching via the superconductor. Than we have the occurrence of non-zero 
vector potential as given by and order parameter of the form 

. As before we fix gauge VA = so we have V x V x A = — V^A(r) = 
-{£^ + 7ir)Mr) = ^mr)?)Mr) and = ^M^r)). We have the GL 
equation as for <^ = 0, a|V(r)| + m{r)? - ^{£^ + l£mr)\+ c^A^{rf = 
Since the system has the rotational symmetry than the solution of the given 
equation is valid for all cases of 4>. 

We can go further and describe the occurrence of the vortex for the 
two dimensional system as by assuming the dependence of the il){r,z,(l)) = 
|V'(r, 2)|e*'^^*('^.^) 

V X V X A = -W'A{r) = + + £,)A^{r) = ^(|^(r)|2)A4r) We 

have a\i,{r,z)\+mr,z)\^ _ + + £,)\i,{r,z)\ + c:,A^{T,zf = 

Having those equations we can describe the properties of the systems as 

finite thickness superconductor with magnetic field perpendicular to its surface. 

We can also describe the vortex in the superconductor-ferromagnet system or 

vortex in the superconductor-normal metal system. The relaxation method can 

be used to solve the given set of differential equations. 

One of the interesting cases to be studied is given below. 
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Similar case to the given above can be accounted by the simple algorithm 
described earlier. What we need to do it to have the square symmetry of the 
system. 
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Another interesting case is given as below 
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3.10 Tunable Josephson junction as the tool for detection 
of the electromagnetic radiation or for other appli- 
cations 

It is possible to produce the tunable Josephson junction using the ferromagnetic 

strip by connecting uJJ ferromagnetic to the magnetically biasing environment. 
If we control externally the state of the magnetization of the ferromagnetic bar 
we might produce the Josephson junction whose sensitivity to the external radi- 
ation can be tuned quite continuously. Also the tunning of the magnetization of 
the ferromagnetic bar on the top of superconductor changes RCSJ parameters 
of the circuit. Therefore we might obtain tunable Josephson junction circuits 
what can have its importance in future THZ electronic circuit design. 

3.11 Vortex role in the detection and emission of the elec- 
tromagnetic radiation 

Basically we expect that the pinned and unpinned vortices are present in the 
system. Because they have the vortex core that has non-superconducting phase 
component therefore they can emit and absorb radiation. Vortex movement 
in the superconductor is the main source of the dissipation and therefore con- 
tributes to the observed electric resistance both in ac and in dc case. 

In the system we might have both Abrikosov and Josephson vortices. As 
pointed by Franco Nori, movement of the Josephson vortices creates more dis- 
sipation than movement of the Abrikosov vortices. 

Also bigger core means bigger cross-section for the absorption of the radia- 
tion and bigger cross-section for the emission of the radiation. The dynamics of 
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vortex movement in uJJ due to the system complexity shall be rather modeled 
phenomenologically. 

Certain simplifications of vortex movement are due to the possible reduction 

of dimensionality as to one or better two dimensions. This is the case of the 
Josephson vortices. However certain effects as switching between Josephson to 
Abrikosov vortices or reversely can be done only in the case of 3 dimensional 
considerations. 

Also pancake vortices are possible to occur in the system and they are a kind 
of non-linear combinations of Abrikosov and Josephson vortices. 

3.12 Presence of single and triplet component of order 

parameter 

It shall be underline that the presence of the ferromagnetic bar shall induce the 
existence of triplet component in the neighborhood of it. In the case of singlet 
order parameter the magnetic field is the factor which destroys the Coopair 
pairs, which are the carriers of the superconductor state. Such situation is not 
occurring in the case of the triplet superconductivity. It might turn out that 
the coming singlet Cooper pair is partly converted to triplet Coopair pair, the 
triplet Coopair pair is then partly converted to singlet Coopair pair as well. In 
such case, we shall observed bigger Josephson current than those expected from 
the s-wave tt Josephson junctions. Please refer to the figure 12. 

3.13 The methodology to compute the Josephson current 
3.13.1 The tunneling Josephson current 

The most known type of Josephson junction is the tunneling Josephson junction 
of rectangular geometry. This is due to the fact that for the simple geometry 
we can describe simple tunneling model. 

We have given the Hamiltonians of the left Hl and right Hji and supercon- 
ductor as 



where a-i are the Pauli matrices. 

In general case we have the Hamiltonian H describing the Josephson junction 
in the form as 




H = Hl + Hfi + Ht 



Here the tunneling Hamiltonian is given as 



Hermicity of Ht gives {T^^y = T^f . Therefore we have 
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k,q 

where the Hamiltonian accounting for the creation(annihilation) of particle 
on the left side superconductor is (d) and on the right side is c' (c). 

The Tj^* is the tunneling matrix between left and right superconductor. 

We might compute the Josephson junction properties by the computation of 
the partition function defined as Z = Tr{e~~). 

In the Hamilton formalism we have the relations between conjugate variables 
q and p given as 

We have the following relation between free energy F and Hamiltonian as 
given 

1 , d d ^ 

-^^^"(^)=<^^>=^^ 

Josephson current can be expressed in the way as 

Wc have the conjugate variable given as the phase of the order parameter (j) 
and the operator giving the numbers of the particles as N and we have [(j), N\ = 
ih, where [] is the anticommutator. 

We can also account for the existence of the spin-current in the system. 

We define Sk = J2q,a,i3 4+q,a^»,pCk,i3 

Here angle 8„ is the angle in the plane perpendicular to the direction n. 
Similar as before we have [6, Sn] = ih 

The Dyson equation is given as = Gq ^ + S and hence G normal prop- 
agator can be written as 

G«Go + GoEGo + GoSGoSGo + ... (2) 

where Go is the normal propagator without presence of the order parameter 
and S is the self-field, which accounts for the presence of the order parameter 
(self-energy) so we have 



Go 



g \ / A 

-9*1 ' I Al' 
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The normal propagator structure becomes more complicated in the case of 
occurrence of the magnetic field as 




) 



The presence of magnetic field can cause the spin-flop during tunneling pro- 
cess. Therefore we shall consider the following structure of the tunneling matrix 
as 



For Sc-I-Sc (Superconductor-Isolator-Superconductor) Josephson junctions 
with no impurities and small tunneling current the tunneling matrix is the 
identity matrix. However for the SIS junctions with some magnetic impurities 
in I region the tunneling matrix is no longer identity. It is not diagonal in the 
case of S-Fe-S Josephson junctions. 

It should be underline that in the proximity of the singlet superconductor 
to the ferromagnetic element the triplet superconducting phase shall be in- 
duced. Therefore we account the region of the superconductor very close to 
ferromagnetic element and also the region of the ferromagnet very close to the 
superconductor. 

Then the form of the self-energy shall be changed and should be as 



Here M accounts for the magnetization of the superconductor and A ele- 
ment consist both singlet and triplet component. The triplet component of the 
superconductor can be characterized as 

A = dto a=iAtay where dt stands for the vector denoting 3 components of 
the triplet phase. In general case we ca have the mixture of single and triplet 
case and express A = (ifO(T + (i5l2x2. Here ds is the complex scalar field. 

We can also express the tunneling Josephson junction in terms of the action 

as 



where Sl,Sl describes the action for the isolated left, right superconductor and 
St describes the action connected with tunneling process. 

For non-interacting fields we have the inverse of normal propagator Gq ^ 
given as 



Here w„ is the Matsubara frequency and Cfe is the kinetic energy of fermion 
gas. 





S — Sl + Sii + St 



Gq ^ = diag{-iujn + efe, -iw„ -|- efe, -ia;„ -|- e_fe, -iuJn + e_fe) 
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Wc introduce fermion coherent states C and performing a Hubbard-Stratonovich 
transformation, we arrive at an effective action (in Euclidean time) , which reads 
(summation over spin indices is given by a,(5) 




+A(a,;3)fc,CV,/L,„- E A(fc',g)l_^l/(fc,fc')"'A(fc,g)0,„ 

k,k\q 

We have the following property 
where -D[C,C^] accounts for the path integral. 

3.13.2 The model tt Josephson junction in the second quantization 
scheme 

Since uJJ have some similarity to the tt Josephson junction wc would like to 
briefly review its physical properties. We are given the system as sci — /e — SC2 
with thin ferromagnetic layer. The presented model is one dimensional. Than 
we have the Hamiltonian of the system given as 

H = Hsci + Hsc2 + Hfg + Hf 

where the Haci,Hac2 and Hf^ are the Hamiltonians of the non-interacting scl, 

sc2 and fe. The interaction is being accounted in the tunneling term Ht that has 
two components H^l and i?t2- We have expression for both values of (j=l,2) as 

k,p,q k,p,q 

Here Tf^f or T^^ represents the tunneling matrix and Cjk,a is the annihilation 
operator of an electron with momentum k and spin a in superconductor j and 
a'p ^ is the creation operator of an electron in fe. 

Let us assume that the ferromagnetic bar is magnetized in parallel to the 
superconductor surface. This implies that the electrons in the Coopair pairs 
close to the fe surface with align in the opposite direction. 

Tunneling of electron with opposite spin is more difficult than for the case 
of the electron with spin of the same orientation. However, if the electron with 
opposite to the fe magnetization spin emits magnon than it can participate in 
the tunneling process much more easily. 

Let us denote the Fourier transform of the localized spin operator with mo- 
mentum q as S*' = {S^, S^,Sl) in the ferromagnetic material. 
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In magnon- assisted tunneling process, a down-spin electron in scl (sc2) pen- 
etrates through the barrier into ferromagnetic part as an evanenscent wave, 
and changes into an up-spin electron by emitting a magnon via the exchange 

interaction. 

The Josephson coupling energy is calculated by the change in the thermo- 
dynamic energy U = -kTlnZ where U = Tr{e-"/^'^). 

Wc can assume that the superconductors 1 and 2 having phases Oi and O2 
are uniform and hence have the anomalous Green functions as 



= Ns{Q)e 



Using the known total coupling energy we compute the Josephson current / 
According to S.Maekawa we have 



where Vm = ^mnkT and xt~iQ) the transverse spin susceptibility with 
momentum q and frequency i/m and (d) is the magnon propagator. 
We have 

^^rw = E'^-"(9K''^"""^ 

is the magnon propagator traveling with a distance d and 

go = (3/2)|rt4,Tt,tP(GT/GK)V(fc/W) 

Here kf is wave- vector value at the Fermi surface, d is the thickness of ferro- 
magnetic material and d is the mean free path for the electrons. 
We have with non-spin flop tunneling 

Gt = {2ne/h)m,^fN,{0)Nf{0) 

and quantum resistance is given as 

e2 

Gk = 

n 

3.13.3 Extension of the tunneling model to the case of the uJJ 

First of all we do not have isolator in the uJJs. Therefore we shall divide it into 
small squares and consider the following Hamiltonian as 

H = Hl + Hr + He + Ht,l-^r,r-^l + Ht,r->^c,c-^r + Ht,r^c,c^r 
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This is how the problem is being described in one dimension. However, we deal 
with at least two dimensions. 

Therefore we need to consider at least 9 cells as given on the picture attached 
below. 



(i-lj+l) 


(iJ+1) 


(i+lJ+1) 


(i-lj) 


(i,j) 


(i+lj) 


(i-l,j-l) 




(i+l,j-l) 



We need to consider the interaction of every cells with another cell. Than 
we have 

We have given the formula for the 

= \hj >< i,j\Es(ij) 

and 

Alternatively we might write the Hamiltonians in the way as 
and 

This model can be easily generalized to the 3 dimensional case. The existence 

of the magnetic field can be accounted by introducing the additional index to 
the creator and annihilator operators and to the transmission matrix as well. 

In a sense the described model accounts for the existence of the Josephson 
junctions in the systems that arc of weak-link type. Unlike in the tunneling 
Josephson junctions here the tunneling energy (tunneling matrix) has value 
comparably to the self-energy. 

Once the Hamiltonian is constructed that it can be represented by two di- 
mensional matrix. The eigenvalues of this matrix are quite easily to be found 
and are the eigenvalues of the sysytem. 

By this method we can account for the possibility of existence of the network 
of the Josephson junctions in the uJJs. 

3.14 RCSJ model 

3.14.1 RCSJ model for the tunneling Josephson junctions 

RCSJ stands for the resistance capacitor shunted Josephson junction. The two- 
fluid model can be used to describe the properties of the tunneling Josephson 
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junction. Phcnomcnologically wc recognize that the existence of two phases: 
superfluid electrons that can move without dissipation and the electrons in the 
normal state. The presence of the last component can be represented as the 
resistance and capacitance connected in parallel. The superfluid component 
can be represented as the perfect Josephson junction, which was introduced by 
Josephson. 

Wc have the following current-phase equation for the Josephson junction 

given as 

ic{t) = iscit) + in{t) = iosinm)) + IM) |e(t) + c^^m 

where Q accounts for the phase difference of the order parameter across the 
Josephson junction. 

3.14.2 RCSJ model for the uJJ 

It should be underlined that there should exist the RCSJ model for the uJJ. 
Having such model we can extend the library of the programs serving for the 
complex circuit design as Spice. In comparison with the RCSJ model for the 
tunneling Josephson junctions, R is nonlinear and strongly depends on the cur- 
rent magnitude and presence of the external magnetic field so we have R(I,B). 

The capacitance can be induced placing the nonzero perpendicular magnetic 
field. 

3.15 Critical temperature of the uJJ 

The presence of the non-superconducting strip on the top of superconductor 

lowers the critical temperature of the superconductor. For thin-layers of Sc-n or 
sc-fe there are some theoretical works describing the dependence of the critical 
temperature on thickness of superconducting-nonsuperconducting layers using 
the Usadel formalism or Eilenberger formalism. 

Let us consider the one dimensional system as the superconductor-ferromagnet 
system of the thickness d of the superconductor and thickness of ferromagnet 
of di. 

We write the Usadel equations for the system. We have very small order 
parameter so we assume that F„ = sm(6„) « 9„ and G„ = cos(0„) « 1. 

Initially I consider only Sc-N system and I assume no presence of the mag- 
netic field. Then we have 



2uinFn{x) 

in the Sc layer and 

2w„F„(a;) 



-D{F^{x)-^Gn{x) - Gn{x)^F„{x)) 
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After the simplifying assumptions we obtain 



2uJnFn{x) = -D{—Fn{x)) + 2A{x) 
and in the N layer we have decoupled F„ functions as 

2uMx) = -D{^Fn{x)) 
We have the self- consistency relation given as 

Here is the critical temperature of the superconductor structure. At outer 
surfaces of the bilayers of Fe (N) or Sc we have 

In the interface between the superconductor and normal state occurring fox 

x=0 we have 

where the (T„ and cFs are the conductivities of the material n and s in the normal 
state. 

We also have another boundary condition given as 

F„(0-) = F,e(0-) - e7^i^n(0) 

Simple algorithm as used to solve two dimensional BdGe equation can be 
used here to solve the Usadel equations in the limit of small order parameter 
size. 

This algorithm is however only limited to the case of thin superconductor 
layer for Sc-N or Sc-Fc structure, which implies that the order parameter present 
in the superconductor is small, what makes the Usadel equation to be linear and 
easy to be solved. For thicker Sc layer such simple strategy of solution is no 
longer possible. 



3.16 uJJ in the external microwave field 

If we consider the uJJ Hamiltonians and when we have the ferromagnetic 
strip in the system we have modification of the magnetic field by H{x, y, z) = 
IJ,Bs{x,y, z) + Bm{t,x,y, z). The last term Bm{t, x,y, z) is sinusoidal time de- 
pendent. Its dependence on the position is rather difficult to be determined. 
Also the phase in the superconductors is time dependent ^-^(f>{t) = V{t) where 
V{t) is of sinusoidal dependence. 
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3.17 SQUID built on the uJJ 



We can built the superconducting SQUID basing on the uJJ. Especially inter- 
esting effects are expected to occur when the two ferromagnetic strips starts 
strongly couple. 
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If the coupling is strong it may cause the occurrence of the topological Meiss- 
ner effect. The external flux of magnetic field cannot enter the interior of the 
SQUID since it would imply the change of the currents flow in the system. This 
process is schematically described below. 
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The similar process shall occur for the SQUID built from tt Josephson junc- 
tions. 

Topological Meissner effect can be induced by the external parallel magnetic 
field applied to the surface of uJJ or tt SQUID. 

The Meissner topological effect (MTE) occurs in 3 dimensions and is rather 
difficult technically to be computed for the GL formalism. However there is 
also possibility of having two dimensional MTE as for two very long sc-fe strips 
connected at the ends of the cross section the same as depicted on the figure 
above. Than investigation of such effect would involve 2 dimensions only as 
Az{x,y), \ip{x,y)\. It can be done with use of 2D GL, BdGe, Usadel or Eilen- 
berger approach. 

The elongated SQUID showing the MTE seen from the top is given below. 
One dimension is assumed to be very big compared to all other dimensions. 
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The structure given below also should exhibit MTE in the case when fe 
element are magnetized around the superconductor. 
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We can also quite easily simulate 2 dimensional SQUID with two uJJ Fe 
or uJJ Ee and Nujj or 2uJJ N or tt— JJ ,uJJ and or single uJJ so on in GL 
formalism. 

3.18 GL equations for uJJ SQUID 

Using certain symmetries of SQUID made of s or d wave superconductors we 
might simply GL equations and use the relaxation method to determine its order 
parameter distribution and vector potential distribution. Let us consider 2 uJJs 
made of non-ferromagnetic strips. 

3.19 Entropy and specific heat of uJJ 

Knowing the entropy wc can quite easily compute the specific heat of the uJJ. 
We can know entropy by solving Usadel equation. 

4 Future perspectives 

The relaxation method seems to be the proper tool to study more complex 
structures as depicted on figure 9 and 10. The obtained solution of the order 
parameter shall be more exact when we plug the order parameter obtained from 
GL solution to the self-consistent solution of BdGe algorithm. 

The suggested geometry of described uJJ shall have the potential to be 
used for THZ Josephson junction. The given geometry shall be theoretically 
investigated in the detail for s, p, d and f wave order parameter in the singlet, 
triple case and the mixed case. What is more, the conduction of the experiments 
is necessary to find the properties of the uJ Js and find their correspondence to 
the existing theories. 

We shall also use Keldysh formalism to model the behavior of uJJ in the 
presence of the external microwave field. It should be underlined that since the 
dimensions of the uJJ available experimentally are about l/um they will give 
the spectrum of energies too closely positioned, so its applications for the qubit 
are rather limited. However one could exploit the existing pinned vortices and 
here the possibilities seems to be more promising. In order to achieve the real 
application of the uJJ to the qubit we need to have the non-superconductive 
strip of width comparable with superconductive coherence length. 
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Fig. 1. The layout of the structure is given below(picture produced by dr. L.Gomez) 
Fig. 2. The scheme of the computed structure of solutions of GL equation for s- 
wave order parameter for 2 dimensional case with the use of relaxation method. 
No currents and no magnetic fields are present in the system. 
Fig. 3. The cross section of the order parameter for the rectangular shape su- 
perconductor of infinite length (upper picture) . The relaxation method used for 
solution GL equation gives the free energy dependence with iterations as given 
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below. 

Fig. 4. The dependence of the total electric current flowing via rectangular 
shape superconductor of infinite height on the iteration step in the GL relax- 
ation method. Please refer to the figure 3. 

Fig. 5. Scheme of 3 basic configurations of magnetic field which are due to the 
certain magnetization of ferromagnetic bar, which is being put on the top of 

superconductor 

Fig. 6. The uJJ magnetic field triangle, which accounts for all possible situation 
of magnetization of the ferromagnetic bar, with given fixed magnitude of mag- 
netic field in geometric center of the ferromagnetic bar. One direction of this 
magnetic field corresponds to one point in this triangle. Changing the direction 
of magnetic field, but keeping its magnitude constant, we move from one to an- 
other point in this triangle. Any point of this triangle is the reference point to 
establish other parameters of unconventional Josephson junction (as resistance 
of uJJ in the limits of small currents or capacitance of the uJJ in the limit of 
small electric current flowing via junction). Points A, B, C correspond to the 
configurations of magnetic field as given on the sub- pictures of Figure 5. 
Fig. 7. The stages of solution after iterations t0,2t0, 3t0, 4tO of the GL equation 
for rectangiflar of infinite height and flnitc a=b dimensions. (please sec below). 
No currents and no magnetic fields are assumed to occur in the system. 
Fig. 8. (Picture below )Distribution of the order parameter for 3D SNS system 
obtained by the solution of GL equation with use of relaxation algorithm. No 
currents and no magnetic field is assumed to occur in the system. Horizontal 
axes correspond to x,y coordinates and vertical axes refers to the magnitude of 
order parameter. z=const for the given picture. 

Fig. 9. The SQUIDs built with use of uJJ. The picture above describes the sys- 
tem showing topological Meissner effect as predicted in the Master of Science 
thesis of Krzysztof Pomorski (2007) at the University of Lodz (Predictions of 
physical phenomena in the mesoscopic superconducting structures and super- 
conductors ). In order to prevent the occurrence of the topological Meissner 
effect one has to decouple the ferromagnetic strips what is depicted on the pic- 
ture below. 

Fig. 10. The possible configuration of the magnetic field in the sc-fe system, 
which exhibits the occurrence of exotic vortex in the superconductor. The stud- 
ies of this structure shall be in the potential of GL relaxation method. Such 
vortices are expected to occur in uJJ system. 

Fig. 11. Order parameter for the superconductor- ferromagnet-superconductor 
system of fixed thickness of ferromagnetic material for different values of mag- 
netic field. We assume the magnetization of the iron to be perpendicular to 
the surface of superconductor. No electric current perpendicular to the super- 
conductor surface occurs in the system. We can observe the oscillations of the 
order parameter in the system as given by Buzdin . This is example of the 1 
dimensional GL problem solved by relaxation method. 

Fig. 12. The behaviour of the Coopair pairs in the proximity of the magnetized 
ferromagnetic bar. One should expect the occurrence of the induced triplet 
phase of superconductor in the proximity of the magnetized ferromagnet. Such 
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situation shall occur in the uJJ as well. 

Fig. 13. Parameterization of the unconventional Josephson junction 
Fig. 14. Possible Spin-flip in uJJ junction. 

Fig. 15. The contour plot of real valued order parameter distribution (no mag- 
netic field and currents occur in the system) in the rectangular of infinite height 
and dimensions of a,b with the rectangular hole inside also infinite height and 
dimensions c,d. The order parameter is obtained from the relaxation method in 
GL formalism. Such system is preludium to study more complex superconduc- 
tor geometries, including 2 and 3 dimensional SQUID built with unconventional 
Josephson junctions. 

Fig. 16. Solution of GL equation for 3 dimensional uJJ with following parame- 
ters:Nx=Ny=Nz=20 (nr of grids elements in 3 dimensions), Lx=Ly=Lz=4 (co- 
herence length=0.5),ybegin=5, yend=15, beginning and ends of normal strip on 
the left side of picture), the number of iterations is 1000. 

Fig. 17. Solution of GL equation for 3 dimensional uJJ with following param- 

eters:Nx=Ny=Nz=20 (nr of grids elements in 3 dimensions), Lx=l, Ly=Lz=4 
(coherence length=0.5),ybegin=5, yend=15, beginning and ends of normal strip 
on the left side of picture), the number of iterations is 1000. 
Fig. 18. The eigenspectrum of 2 dimensional solution of BdGe equations for uJJ 
with no currents and fields. Parameters from Fig. 16 are being taken. 
Fig. 19. The scheme of pi- Josephson junction with long Fe channel. Conse- 
quently one dimensional Eilenberger equation can be solved and it gives non- 
sinusoidal current phase relation. 

Fig. 20. Multiple BdGe wavepacket reflections in SNS system. The red arrows 
stands for hole wavepacket, while black arrow stands for the electron wavepacket. 
Fig. 21. Non-uniform BdGe wavepacket reflection in uJJ. 
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